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A LIFTING  SURFACE  OF  AVERAGE  ASPECT  RATIO  IN  A SUBSONIC  GAS  FLOW 


A.  N.  Panchenkov 


(Kiev) 


For  nany  technical  applications  there  mst  be  information  on  the 
hydrodynaaic  characteristics  of  lifting  surfaces  of  average  aspect 
ratio  of  arbitrary  planforn.  However,  in  developing  an  analytical 
theory  of  a lifting  surface  of  average  aspect  ratio  there  are 
difficulties  involving  peculiarities  of  the  general  problem  of  the 
notion  of  a lifting  surface  in  a gas  flow. 


Methods  of  the  potential  of  accelerations  have  aade  it  possible 
to  obtain  a two-diaensional  singular  integral  equation  for  an 
arbitrary  lifting  surface  under  various  conditions  [1],  Further  study 
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proceeds  along  two  directions. 

I 

1.  Development  of  numerical  aethods  of  solving  the  integral 

equation.  The  most  widely  developed  aethod  in  various  studies 

is  the  aethod  of  collocations,  which  in  this  problem  has  a simple 
hydrodynamic  interpretation.  Distributing  discrete  vortices  along  the 
lifting  surface,  let  us  trace  the  foraalistics  of  the  method  of 
collocations;  here  the  number  of  vortices  determines  the  precision  of 
the  solution. 

k system  of  one-dimensional  integral  equations  of  the 
collocation  method  can  be  obtained,  avoiding  operations  of  the 
methods  of  acceleration  potential  - by  exaaining  the  physical  picture 
of  the  flow,  caused  by  f? -shaped  vortices.  This  method  is  used  in  [ 2 ] 
to  thoroughly  stui y a broad  class  of  stationary  problems  of  a lifting 
surface. 

2.  Development  of  analytical  methods  for  lifting  surfaces  of 

limited  planform.  The  aost  widely  developed  theory  is  that  of 

the  Prandtl  lifting  line  - for  wings  of  large  aspect  ratio  - and  the 
theory  of  low-aspect  wings  (the  Jones  theory)  . 

The  Prandtl  theory  is  "incorrect*1  in  the  limiting  sense  [7];  in 
this  connection  it  is  not  possible  to  construct,  using  familiar 
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methods,  a theory  of  a wing  of  average  aspect  ratio  in  which  the 
Prandtl  problem  would  be  the  zero  approximation.  This  "incorrectness" 
in  the  limiting  sense  of  the  Prandtl  problem  creates  mathematical 
difficulties  and  has  attracted  the  attention  of  researchers.  For 
example,  H.  Van  Dyke  [8]  again  turns  to  this  problem  and  develops  a 
solution  based  on  various  expansions  and  on  the  connection  of  the 
solutions  in  the  problem  with  "singular  perturbations"  (using  his 
terminology).  The  theory  of  a low-aspect  wing  is  "correct"  in  the 
limiting  sense,  and  it  can  be  obtained  by  passage  to  the  limit  from 
the  general  problem  of  a lifting  surface  [7]. 


The  theory  of  a low-aspect  wing  gives  satisfactory  results  only 
for  X « 1,  and  its  extension  to  an  aspect  ratio  of  ~1  and  more  leads 
to  high  errors,  particularly  in  determining  the  moment  on  the  wing. 
However,  its  "correctness"  opens  the  way  for  development  of  methods 
in  which  the  Jones  theory  will  be  the  zero  approximation. 

The  authors  of  [4,5]  derive  a solution  to  the  problem  in  the 
form  of  a series  in  powers  of  X2,  which  does  not  appreciably  expand 
the  boundaries  of  applicability  of  the  theory. 


In  [7]  we  propose  seeking  the  solution  to  ths  problem  in  the 
form  of  a series  in  posers  of  .=  4- 1 - -j  • shich  can  give 

converging  results  for  a broad  range  of  X.‘  The  results  obtained  in 
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this  work  ace  encouraging;  however,  the  unwieldiness  of  the 
asyaptotic  Methods  of  the  perturbation  theory  cbviously  Leads  to 
difficulties  in  daveloping  a theory  for  wings  of  arbitrary  fora. 

In  the  present  work  we  develop  an  analytical  theory  of  a lifting 
surface  with  average  aspect  ratio,  based  on  the  use  of  integral 
operators  [3]  that  convert  the  solutions  of  the  two-diaens iona 1 
Laplace  equation  into  solutions  of  elliptical  equations.  The  general 
representation  obtained  for  the  velocity  potential  includes  all  the 
results  known  froa  the  literature  as  the  zero  approxiaation  (a  wing 
in  plane-parallel  flow,  the  Prandtl  theory,  the  Jones  theory,  and  the 
Laidlaw  theory  of  wings  of  arbitrary  aspect  ratio  [6]).  The  first 
approxiaat ion,  examined  in  this  work,  leads  to  a linear  second-order 
ordinary  differential  equation  for  load  distribution  over  the  wing 
chord.  This  differential  equation  is  easily  solved  for  arbitrary  wing 
planf ora. 

As  applications  of  the  general  results,  in  this  work  we  study 
wings  of  rectangular  and  triangular  planfora.  The  vast  aaount  of 
calculation  data  given  in  [2]  is  in  good  agreeaent  with  the  obtained 
results  for  the  intervals  of  practical  interest  up  to  \ - 4. 

With  an  increase  in  Each  nuaber  the  liaits  of  applicability  of 
the  results  increase,  and  when  H > 0.75  they  are  valid  to  x - 10. 


i 


I 
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This  fact,  in  particular,  leads  to  the  conclusion  that  for  high  Mach 
nuabers  eren  high-aspect  wings  Bust  be  studied  using  the  theory  of 
wings  of  average  aspect  ratio. 

The  problea  of  the  stationary  aoveaent  of  a lifting  surface  in 

an  incoapressible  fluid  reduces  to  the  boundary-value  problea  for  the 
Laplace  equation: 

Aq>  «■«  0; 

?,  - F (*);  g £ «il  (1) 

«px  «*  0;  x.  V.  » £ J„ 

where  s,  is  a surface  aoving  in  a fluid  without  perturbations, 
coinciding  with  the  projection  of  the  lifting  surface  onto  plane  Oyr; 

s,  is  a seai-inf inite  surface  aoving  in  a fluid  without 
perturbations  and  beginning  at  the  trailing  edge  of  surface  st- 


nsing  the  basic  idea  of  Bergaan's  aefchod  [3]  of  the  possibility 
of  representing  the  solution  to  an  elliptical  equation  in  teras  of 
the  solution  to  the  Laplace  equation,  for  solving  the  Laplace 
equation  in  three- diaensional  space  we  can  write  the  representation 


+A.W 


Tfli  (*•  x)  l~n9  (f) 


■,  w <*•  *> + s. 

(V.  *>  + ^ **„  (V.  *)0| 


,(») 


(2) 
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Her*  Vo.  and  *0.  are  solutions  of  the  two-dimensional  Laplace 
eq uation. 

Series  (2)  contains  faniliar  theories  as  lower-order 
approximations: 

3)  * = ♦o  (*#*)  - plane-parallel  flow; 

2)  ♦ - ?(?)♦#  (x,z)  ♦ * ~ Prandtl  theory  of  a high-aspect 

wing; 


3)  P(*Mo(7»*)  “ theory  of  a low-aspect  wing; 

~ F(y)  *0  (x#y)  ♦ P(x)Xg(z#y)  — Laidlaw  theory  of  a wing  of 
arbitrary  aspect  ratio  [6]. 

Let  us  state  the  problem  about  the  realization  of  the  algorithm 
used  to  calculate  series  (2). 

In  coordinates  *,  y. « the  Laplace  equation  will  have  the  fora 


*»  + *-.;  + x' -°- 


(3) 
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The  connection  along  the  coordinates  is  given  by  the 
relationships 

x ■=  xb  (y); 

V - W <*); 

» — 7i>  (y), 

vher*  X (.«)  — is  the  relative  span  in  section  x; 

b(y)  is  the  half-chord  in  section  y; 

li*)  is  the  sesispan  in  section  x. 

If  ths  integral  operator  realizes  the  napping  n't*,  y.  *)  -*  q> (<r.  v.  ’«)  and 
nst  is  ths  solution  to  the  equation 

(4) 

such  an  operator  will  be  called  the  integral  Fourier  operator. 

For  aquation  (4)  there  exists  an  integral  operator  of  the  first 
kind  with  the  degenerate  kernel 

v"»,  + iMy)(v.«»i)^  (5) 
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•her.  ».-c.  Is  the  associated  solution  to  the  lapi.ce  equation. 

"*  easily  find  that  function.  *.  and  f.  ..ti.f,  the  functional 
relation  sh  ips 


- B*Kn  (y)  r„  (v,)  + < (y)  Fn  (yt)  - Kn  (y)  F ' (yt)  ==  0. 


(6) 


V 2[/fn  (y)  <*)]•- Q. 

n«*  I 


(7) 


Let  operator  T be  given  in  the  fo 


ra 


rV~2 k!  V (*’  y-  7)  N <*' 


(8) 


then  froa  (3)  ve  g*t 


V h V-.  I-  - **'  1*  Uf  _ n.  <J  , . N’ 

W **  ' /V  (ft,  x) * H - °'  *■  M -fir 


- ~ *’  (*)■ 


If  function  M(h.r)  is  the  solution  to  the  differential  equation 


A"  (ft,  x)  X*  (i)  + fl*  (*)  \ (ft,  X)  = o. 


(8) 


operator  T in  fora  (8)  realizes  the 


required  napping. 
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Nov  ve  have 

q>  = TV  - ji-  j '•'o  (A,  y,  i)N(k.x) „ . -f 
0 

00  ^ ^ 

+ 2F  2 j J *«  w v.  <*•  **>  * <*•  7> (y“>  diTt-  ( 1 0) 

n**l  ® 

Pron  fora  (10)  ve  caa  obtain  the  second  part  of  series  (2),  if 
we  set 

oo 

2~  j % (*.  y,.  *)  'n  (*.  2)  <<*  = pm  w *>*  6 *7- 

The  first  part  of  series  (2)  can  be  obtained  in  the  sane  Banner, 
introducing  the  operator  Ta  that  realizes  the  napping 

V (x,  k,  z)  - 9 (x,  y,  i). 

Nov  let  us  apply  forn  (10)  to  the  problen  of  the  motion  of  a 
lifting  surface. 

If  we  retain  only  the  first  tern  of  series  (10),  we  get  an 
approximation  corresponding  to  the  theory  of  a low-aspect  wing.  I 

In  accordance  with  this,  ve  assume 

» +i  <?r  (x.  ii) 

srj  r‘u-  * - srrfej 

+1  _ _ . 

where  r (x,  n)  = J vd-  n)4  i*  circulation  about  the  contour  L € H + *)  in 


(1!) 


w 
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auction  i), 


Let  us  introduce  in  (11)  an  assumption  typical  for 
oaly  of  the  class 

r (*•  i)  =*=  r (x)  y (n).  v (0)  = i, 


then 


+i  - 


a problem 


"V*.  y.  0)  = f J (nl-dnG  (k)- 

J,  to  - n) 


(12) 


Formula  (11,  defines  the  r.lme.  of  the  induced  velocity  along 
axis  z on  the  lifting  surface.  Defining  the  induced  velocity  from 

(10),  we  get  the  one-dimensional  singular  integral  equation  of  the 
problem 

-L  r,Mrr«  1 i v r f*(Ft)  ? „ - 

2*  J v (n,U w + -e, J (-^j i *« <*  *> ° w * (*. *) x 


X d/My, 


dri  = f (x.  yj. 


(13) 


Let  us  e,a.i„e  the  pcobl..  1.  th.  necond  approxin.tion. 

retaining  i„  (13)  only  the  terns  containing  r,  * 1 ,.d  x,  ® _ ‘ 0.wl, 
Prom  (9)  we  have 


and  the  equation  will  have  the  fora 


<15> 

Equation  (15)  can  be  solved  by  various  approximate  methods.  In 
particular,  this  equationis  regulated  by  the  action  of  the  inversion 
operator  in  the  class  of  functions  unbounded  at  the  ends,  and  for  its 
solution  we  can  use  methods  of  solving  the  Fredholm  equations.  Let  us 
give  the  solution  to  equation  (15)  for  a lifting  surface  with  a 
constant  angle  of  attack  over  the  span  in  the  variational 
approximation,  assuming  y (n)  - K 1 — V.  Operating  on  both  sides  of  the 

HI  

equation  with  the  operator  1- f Kt  -t-y*dy  , we  get  the  differential 

1 -i 

equation 


r (*)  i1  (*) 
Tlx)  8 


F(x). 


(16) 


Equation  (16)  is  easily  solved  with  arbitrary  functions  X(x)  and 


F (x) . Thus  we  can  study  a broad  class  of  wings  with  varying 


I 
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planf  orms. 

As  the  first  example  of  using  this  developed  theory,  let  us 
examine  the  problem  of  a plane  lifting  surface  of  rectangular 
planform,  for  which  we  have  very  detailed  results  obtained  by  the 
collocation  method  in  [2], 

In  this  problem  equation  (16)  will  have  the  form 

r (*)  — 2?  r (x)  = x.  ( 1 7) 

8 

The  solution  to  (17) , containing  two  unknown  constants,  will  be 

s YT  _ 3 vr 

r(*)  = x[i  + d|«  ' + iV  x x.  • 

As  one  of  the  conditions  for  determining  the  constants,  let  us 
take  the  condition  of  the  Joukowski'Chaplygin  postulate: 

p (_  |)  » o.  P 8) 

The  second  condition  can  be  the  condition  of  the  boundedness  of 
the  solution  with  i -*  00  , or  the  condition 


r (i)  - o. 


(19) 
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vr 


Pro*  condition  (It)  A% . T end  the  solution  is  transfor»ed  to 
the  fora 


rw 


r _Lir  g/?  i 

— X 1^1  — ch  — X — (X  4-  1)1. 


(20) 


Proa  the  condition  of  boandedness  of  tho  solution  when  * 
while  froi  condition  (18) 


« » >r 

— « t 


» »T 


ch 


i V 2 


Now  let  is  define  the  conditions  on  the  wing. 

The  lift  and  nonent  coefficients  of  the  wing  are  defined  by  the 
fornulas 


c»  “ -y-  ar  (—  I). 


r „ 

Cm  - -jr  “ 


| frwrf,- re- 


using condition  (18) , we  get 


(21) 


(2  2) 


(2  3) 
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e _™Lx 
* If 


! + • 


ch 


a yf  2 kit 


(24) 


The  position  of  the  aerodynamic  center  of  the  wing  is  defined  by 
the  formula 


2 


2 Vi 


ch 


A V2 


*f-  1 


I — 


ch 


0^2 


(25) 


If  we  determine  A,  from  the  condition  at  infinity,  the  lift 
coefficient  will  be 


C,-T[l 


— T j , 


(26) 


Now  we  can  compare  the  results  obtained  with  the  familiar 
calculation  data. 

In  Fig.  1,  curvs  1 represents  the  results  obtained  by  S.  n. 
Belotserko  vskiy  for  £ 2 Js  curve  2 corresponds  to  calculation  by 
formula  (2  6);  curve  3 is  for  calculation  per  formula  (23). 

Figure  2 shows  a graph  of  *r  vs  A.  curve  1 represents 
calculation  per  formula  (25);  curvs  2 - the  results  of 
Be  lot  serkovskiy 's  calculations  [ 2 J.  As  Fig.  1 shows,  formulas  (22) 

and  (23)  agree  well  with  data  for  calculation  uf  to  aspect  ratios  of 

~ V. 
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He  know  that  the  theory  of  a low-aspect  wing  gives,  for  a 
rectangular  wing  (see  Fig.  1),  coincidences  with  the  precise  results 
up  to  X - 1 with  respect  to  lifting  force,  but  leads  to  great  errors 
when  determining  the  moment,  even  for  X - 0.2.  Therefore,  conformity 
of  the  developed  theory  with  calculation  data,  with  respect  to 
moment,  is  of  great  interest.  From  Fig.  2 we  see  that  even  formula 
(25)  reproduces  well  the  displacement  of  the  aerodynamic  center  of 
the  wing  in  the  set  X£(o  + <|  . Formula  (26)  agrees  well  with  the  data 
in  monograph  [2]  for  x £ |o  + «|  . However,  determination  of  Aa  based  on 
the  condition  at  infinity  leads  to  great  errors  in  determining  the 
moment. 


For  the  subsonic  flow  of  a compressible  fluid,  in  fornulas 
(20)  — ( 26)  it  is  necessary,  in  place  of  X,  to  set  x***xp  and  P - I i ~ AT* , 
except  for  the  X appearing  with  the  multiplier  in  (21)-(26). 

It  is  interesting  to  note  that  with  increasing  Nach  number  the 
set  in  which  the  obtained  results  are  valid  expands,  and  already  when 
H - 0.75  xe  |6  + 10)  . 

Now  let  us  examine  the  problem  of  a flat  delta  wing.  The 


function  X (x)  for  a delta  wing  has  the  form 
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M*)  - y |*| 

and  the  general  solution  to  equation  (16)  will  be 


" 2 1 °* 2 • (2  7) 

We  deteraine  the  constant*  B,  fro«  the  conditions 


(SLr0' 

(28) 

mi.- 

(29) 

n at2~a‘ 

0,2*'  +a,2(|-a>>  ' fl*“— C + 2fl|)- 

(30) 

The 

foraula 

lift  coefficients  for  a delta  wing  will  be  defined 

by  the 

where 

* 2 

X®  - 2X. 

T - I + fl, 2*'  + fl,2— .. 


(31) 
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Por  the  aoaent  coefficient  the  following  formula  is  valid: 


_ n 


W 

J r (JT)  </jc  — 2r  ( — 
-2 


(32) 


Calculating,  we  get 


c ® 


l--AT2‘,'+s+il.2->. 

“i  + 2 a. 


(33) 


(34) 


The  location  of  the  aerodynaaic  center  of  the  wing  is  deterained 
by  the  foraula 


v, 


(35) 


In  Pig.  3,  curve  1 representes  the  results  of  Lawrence's 
calculations  [ 0 ],  while  curve  2 shows  the  calculations  per  foraula 
(31). 

figure  4 shows  the  curve  *f  calculated  per  foreala  (35).  The 
data  obtained  for  a delta  viag  give  values  for  K*  < 2.5  that  are 
satisfactory  for  technical  applications. 


’-•'V-  - -,v  ' v- 
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The  influence  of  the  Mach  number  in  the  delta-wing  problem  is 
calculated  in  the  same  manner  as  in  the  previous  problem.  Using  the 
developed  method  we  can  also  study  the  problem  of  a lifting  surface 
in  a supersonic  flow. 
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